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Abstract 

In a previous publication, we have established a collinearly-improved version of the Balitsky-Kovchegov 
(BK) equation, which resums to all orders the radiative corrections enhanced by large double transverse 
logarithms. Here, we study the relevance of this equation as a tool for phenomenology, by confronting 
it to the HERA data. To that aim, we first improve the perturbative accuracy of our resummation, by 
including two classes of single-logarithmic corrections: those generated by the first non-singular terms in the 
DGLAP splitting functions and those expressing the one-loop running of the QCD coupling. The equation 
thus obtained includes all the next-to-leading order corrections to the BK equation which are enhanced by 
(single or double) collinear logarithms. We then use numerical solutions to this equation to fit the HERA 
data for the electron-proton reduced cross-section at small Bjorken x. We obtain good quality fits for 
physically acceptable initial conditions. Our best fit, which shows a good stability up to virtualities as large 
as = 400 GeV^ for the exchanged photon, uses as an initial condition the running-coupling version of the 
McLerran-Venugopalan model, with the QGD coupling running according to the smallest dipole prescription. 
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1. Introduction 

The wealth of data on electron-proton deep inelastic scattering collected by the experiments at HERA 
over 15 years of operation has allowed for stringent tests of our understanding of high-energy scattering 
from first principles. This refers in particular to the ‘small-a;’ regime where perturbative QCD predicts a 
rapid growth of the gluon density with increasing energy (or decreasing Bjorken x), leading to non-linear 
phenomena like multiple scattering and gluon saturation [DE]. The simplicity of the dipole factorization 
for deep inelastic scattering at high energy Oil] has favored the emergence of relatively simple ‘dipole 
models’, in which the high-density effects are efficiently implemented as unitarity corrections to the cross- 
section for the scattering between a quark-antiquark dipole and the proton. Such models allowed for rather 
successful fits to the small-a: HERA data at a time where the theory of the non-linear evolution in QCD was 
insufficiently developed and the pertinence of gluon saturation for the phenomenology was far from being 
widely accepted. The first such model — the “GBW saturation model” laisi — provided a rather good 
description of the early HERA data for the inclusive and diffractive structure functions at a; < 10“^ with 
only 3 free parameters. This success inspired new ways to look at the HERA data, which in particular 
led to the identification of geometric scaling |D. The subsequent understanding iHHini of this scaling from 
the non-linear evolution equations in QCD — the Balitsky-JIMWLK hierarchy [iiHni and its mean field 
approximation known as the Balitsky-Kovchegov (BK) equation [18j — has greatly increased our confidence 
in the validity of the pQCD approach to gluon saturation as a valuable tool for phenomenology. 
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Over the next years, new ‘dipole models’, of increasing sophistication, have emerged. On one hand, they 
were better rooted in perturbative QCD, thus reflecting the overall progress of the theory [T9H23]. On the 
other hand, they were better constrained by the advent of new data at HERA, of higher precision. Finally, 
they extended the scope of the ‘saturation models’ to other observables, like diffractive structure functions 
and particle production in heavy ion collisions. Such extensions required more elaborated versions of the 
dipole model, including impact parameter dependence and heavy quarks [13 US] . 

For several years, the theory of high-energy scattering with high gluon density was known only to leading 
logarithmic accuracy in pQCD, which is insufficient for direct applications to phenomenology. For instance, 
the essential running coupling corrections enter the high energy evolution only at next-to-leading order (see 
below for details). To cope with that, the first generations of dipole models involved phenomenological 
parametrizations for the dipole amplitude, which were rather ad-hoc, albeit sometimes inspired by solutions 
to the BK-JIMWLK equations. For instance, the ‘IIM’ fit in m attempted to capture some general features 
of the non-linear evolution, like geometric scaling with an anomalous dimension and the BFKL diffusion, 
that were expected to hold beyond leading order ISIES]. However, the situation has changed in the recent 
years, when the next-to-leading corrections to the BK and JIMWLK equations have progressively become 
available [501 - 155] . This opened the possibility for new fits in which the evolution of the dipole amplitude 
with increasing energy is completely fixed by the theory and only the initial condition at low energy still 
requires some modeling involving free parameters. In that respect, the situation of modern ‘dipole fits’ 
becomes comparable in spirit to that of the more traditional fits based on the DGLAP equation. 

So far, this strategy has been applied [HHHISS] only at the level of the “running coupling BK equation” 
(rcBK) — an improved version of the LO BK equation which resums all-order corrections associated with 
the running of coupling, with some scheme dependence though |30H32j . These corrections are numerically 
large, since enhanced by a large transverse (or ‘collinear’) logarithm, and their resummation within the BK 
equation has important consequences on the evolution — it significantly slows down the growth of the dipole 
amplitude with increasing energy [3IM1I37]. This last feature was indeed essential for the success of the 
HERA fits based on rcBK [5TH551 [55] . The state of the art in that sense is the “AAMQS” fit in [^ , which 
provides a good description of the most recent HERA data [55] (the combined analysis by HI and ZEUS 
for the ep reduced cross-section, which is characterized by very small error bars), with a number of free 
parameters which varies from 4 to 7 (depending upon whether heavy quarks are included in the fit, or not). 

However, the running of the QCD coupling is not the only source of large (but formally higher-order) 
perturbative corrections to the LO BK, or JIMWLK, equations. Besides the running coupling corrections, 
the full next-to-leading order (NLO) corrections to the BK equation, as computed in [55], feature other 
contributions which are enhanced by potentially large, single or double, transverse logarithms. Such terms 
were indeed expected, given our experience with the NLO version [59ll44] of the BFKL equation [45II47] 
(the linearized version of the BK equation valid when the scattering is weak). The NLO BFKL corrections 
are numerically large and thus render the small-a: evolution, at LO and NLO, void of any predictive power. 
There is no reason to expect this problem to be cured, or even alleviated, by the inclusion of the non-linear 
terms describing unitarity corrections |291148] : the collinear logarithms are generated by integrating over 
regions in phase-space where the dipole size is small and the scattering is weak. This has been indeed 
confirmed by the first numerical study of the NLO BK equation [49], which showed that the evolution is 
unstable (the scattering amplitude decreases with increasing energy and can even turn negative) and that 
the main source for such an instability is the large double-logarithmic correction. 

This difficulty reflects the existence of large radiative corrections of higher orders in which formally lie 
outside the scope of the high-energy evolution (since generated by the transverse phase-space), but in practice 
spoil the convergence of the perturbation theory and hence must be kept under control via appropriate 
resummations. In a previous publication [50], we have devised a resummation scheme which deals with the 
largest such corrections — those where each power of Ug is accompanied by a double transverse logarithm. 
Our strategy relies on explicit calculations of Feynman graphs and results in an effective evolution equation 
— a collinearly improved version of LO BK equation —, in which both the kernel and the initial condition 
receive double-logarithmic corrections to all orders. This scheme differs from the ‘collinear resummations’ 
previously proposed in the context of NLO BFKL [5TH55] in that it is explicitly formulated in the transverse 
coordinate space, rather than in Mellin space, and hence it is consistent with the non-linear structure of the 


2 


BK equation. Besides, our equation is local in ‘rapidity’ (the logarithm of the energy, which plays the role 
of the evolution variable), a property which in this context is rather remarkable since the physics behind 
the double collinear logarithms is the time-ordering of subsequent, soft, gluon emissions, which is genuinely 
non-localQ The first numerical studies of this collinearly-improved BK equation demonstrate the essential 
role played by the resummation in both stabilizing and slowing down the evolution 

In this paper, we shall provide the first phenomenological test of our resummation scheme, by using it in 
fits to the inclusive HERA data. To that aim, it will be important to first extend this scheme to also include 
the single transverse logarithms which appear in the NLO correction to the BK equation — that is, the 
NLO terms expressing the first correction to the DGLAP splitting kernel beyond the small-a; approximation 
and those associated with the one-loop running of the coupling. Indeed, such single-log effects must be kept 
under control to ensure a good convergence of the perturbative expansion. Besides, the inclusion of running 
coupling effects is essential for the description of the data, as well known. 

The resummation of the DGLAP logarithms to the order of interest turns out to be rather straight¬ 
forward: it amounts to adding an anomalous dimension (a piece of the leading-order DGLAP anomalous 
dimension) to both the resummed kernel and the resummed initial condition. For the running coupling 
corrections, the situation turns out to be more subtle since, strictly speaking, they cannot be encoded into 
an equation which is local in rapidity. This being said, and following the standard strategy in the literature, 
we shall propose various schemes for introducing a running coupling directly in the local evolution equation 
and test these schemes via fits to the HERA data. 

After these additional resummations, we are led to a new, more refined, version for the ‘collinearly im¬ 
proved BK equation’, namely Eq. below, which will be our main tool for phenomenology. By construction, 
this equation resums the double-logarithmic corrections completely — meaning to all orders in dg = asNc/n 
(ttg is the QGD coupling and is the number of colors) and with the right symmetry factors —, whereas 
the single-logarithmic terms are resummed only partially (but in such a way to include the respective terms 
to NLO). It is rather straightforward to extend our resummed equation to full NLO accuracy, by adding the 
remaining corrections of 0{a^), as computed in [33]. But the ensuing equation would be very cumbersome 
to use in practice, due to the intricate, non-local and non-linear, structure of the pure corrections. In 
this first analysis, we shall adopt the viewpoint that the most important higher-order contributions (say, in 
view of phenomenology) are those enhanced by collinear logs, as explicitly resummed in Eq. ([^, and that 
the pure aj. effects are truly small and can be effectively taken care of via the fitting procedure. A similar 
viewpoint has been advocated in previous fits based on rcBK, but given the importance of the collinear 
logarithms, this assumption was not so well motivated and led indeed to some tensions in the respective fits, 
as we shall later explain. 

Using numerical solutions to this collinearly improved BK equation together with suitable forms for the 
initial condition, we have performed fits to the HERA data for the ep reduced cross-section |38| at x < 10“^ 
and where the upper limit virtuality of the exchanged photon is varied between 

50 GeV^ (a common choice in small-x fits) and 400 GeV^. These fits show several remarkable characteristics. 

(i) The fits are indeed successful: for = 50 GeV^ and two types of initial conditions — GBW-like 
[5| and the running-coupling version of the McLerran-Venugopalan (rcMV) model [SSj — , we obtain a 
per number of data points around 1.2 with only 4 free parameters. 

(ii) The fits are also very discriminatory: they clearly favor some initial conditions over some others, 
and some prescriptions for the running of the coupling over the others. For instance, the standard MV initial 
condition, which truly corresponds to a fixed coupling, appears to be disfavored, whereas a more physical 
version of it, including a running coupling, works quite well. The latter works also better than the GBW 
initial condition, in the sense that it provides a fit which remains stable up to = 400 GeV^. 

(iii) Our fits alleviate some tensions (in terms of physical interpretation) which were visible in previous 
first based on rcBK [HHisiiss] and could be attributed to the choice to replace all the NLO corrections 
with the running of the coupling alone (see also the related discussion in |13|). Notably, our fits prefer 
prescriptions where the QGD coupling as{pi^) is running according to the smallest dipole size, they do not 


^In fact, a non-local equation to resum the double logarithms has been proposed in m- 
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require any artificial ‘anomalous dimension’ in the initial condition, and treat the heavy quarks on the same 
footing as the light ones, in agreement with general expectations from the dipole factorization. 


2. The NLO BK equation and large transverse logarithms 


To motivate the resummations that we shall later perform, let us first explicitly exhibit the large trans¬ 
verse logarithms which appear when computing the NLO corrections to the BK equation [5TH55] . We recall 
that the BK equation describes the rapidity evolution of the S'-matrix S^y = 1 — T^y for the scattering 
of a color dipole with transverse coordinates {x,y) off a hadronic target. The dipole scattering amplitude 
Txy is small in the regime where the target is dilute, but it approaches the unitarity (or ‘black disk’) limit 
Txy — 1 when the target is dense. The separation between these two regimes is controlled by the saturation 
momentum Qs(Y)-, which increases with the rapidity difference Y between the projectile and the target. 

Neglecting the terms suppressed in the limit of a large number of colors TVc 3> 1, one finds a closed 
equation for Sxy, whose NLO version reads as follows [55] 
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where Nf is the number of flavors, b = (lllVc ~ 2fVf)/12A^c; and dg = asNc/n, with the QCD coupling ag 
evaluated at the renormalization scale 

There are two main changes in the structure of the evolution equation as we go from LO to NLO. 
First, the term with a single integration (SI) over the transverse coordinate z only receives a correction of 
order 0{dY) to the kernel, which in particular contains the running coupling corrections proportional to b. 
Second, there are new terms, of order 0(11^), which involve a double integration (DI) over the transverse 
coordinates u and 2 : and which refer to partonic fluctuations involving two additional partons (besides the 
original quark and antiquark) at the time of scattering. The first such a term, which is independent of 
iVf, represents fluctuations where both daughter partons are gluons. The S'-matrix structure therein, that 
is, SxuSuzSzy — SxuSuy, corresponds to the following sequence of emissions: the original dipole (x,y) 
emits a gluon at m, thus effectively splitting into two dipoles {x,u) and (u,y); then, the dipole {u,y) emits 
a gluon at 2 :, thus giving rise to the dipoles {u,z) and {z,y). The ‘real’ term SxuSuzSzy describes the 
situation where both daughter gluons interact with the target. The ‘virtual’ term —SxuSuy describes the 
case where the gluon at 2 ; has been emitted and reabsorbed either before, or after, the scattering. This 
negative ‘virtual’ term subtracts the equal-point contribution [z = u) from the ‘real’ piece, ensuring that 
the potential ‘ultraviolet’ singularity associated with the factor 1 /(m — 2 ;)'^ in the kernel is truly harmless. A 
similar discussion applies to the second DI term, proportional to A^f, except for the fact that the additional 
partons at the time of scattering are a quark and an antiquark. 

In principle, one should be able to undertake the task of solving the NLO BK equation. The hope 
would be that the solution would only add a relatively small correction to the LO result. However, this is 
not the case since there are terms in the kernels of the NLO equation which can become large in certain 
kinematic regimes and thus invalidate the strict ds-expansion. One obvious class of such terms contains the 
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corrections proportional to b in the SI term in Eq. Q, which by themselves bring no serious difficulties: as 
well known, these corrections can be absorbed into a redefinition of the scale for the running of the coupling, 
which thus becomes a dynamical scale (see Sect. below for details). Here, we would like to focus on the 
corrections enhanced by ‘collinear logarithms’, that is, logarithms associated with the large separation in 
transverse sizes (or momenta) between successive emissions. These corrections become large only in the 
weak-scattering regime where all the dipoles are small compared to the saturation scale l/Qs{Y) and the 
equation can be linearized w.r.t. to the (small) scattering amplitude T. This in particular means that one 
can ignore the last term, proportional to Nf/Nc, in Eq. Q since this term vanishes after linearization, as 
one can easily check (by also using the symmetry of the kernel under the interchange u ^ z). 

To be more precise, let us consider the strongly ordered regime 

1/Qs » I 2 : - a;| ~ |z - y| ~ |z - mI » |m - a;| ~ |m - y| » |a; - y\, (2) 


that is, the parent dipole is the smallest one, a gluon is emitted far away at it, a second one even further at 
z, but with all possible dipole sizes remaining smaller than the inverse saturation momentum. Whenever 
appropriate, we will denote by r, u and z the size of the parent dipole, the size of the dipoles involving u 
and the size of the dipoles involving 2 :, respectively, with z^. By inspection of the SI piece in the 

NLO BK equation, it is quite obvious that the dominant NLO term is the one involving a double transverse 
logarithm (DTL), that is, the last term within the square brackets. Still within this regime ([^, we can 
approximate the scattering matrices in the SI term as follows: S^zSzy — — —Txz — Tzy + T^y — —2r(2), 

where the second approximate equality follows since the dipole amplitude for a small dipole is roughly 
proportional to the dipole size squared. Notice that the net result in the approximation of interest fully 
comes from the ‘real’ term, which involves the large daughter dipoles. 

What is not immediately obvious is the presence of a single transverse logarithm (STL) coming from the 
DI term. Let us isolate here the relevant part of the kernel. 
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To implement the limit in Eq. ([^ we can successively write the expression in Eq. ([^ as 
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with (j) the angle between r and any of the two dipoles involving u. To obtain Q, we have first set all dipole 
sizes which include 2 equal to each other, since any subleading term would be suppressed by inverse powers 
of 2 . Then the only 2 dependence left is the one explicit in the prefactor. We have subsequently taken 
the limit r <gC u (by expanding the logarithm to cubic order) and we have finally averaged over the angle 
(j) between the parent dipole and those involving u. Notice that the would-be leading term, of order 1 / 2 ^, 
has cancelled out in these manipulations. The first non-vanishing term, as visible in the r.h.s. of Eq. Q, 
is suppressed by r'^/v?, thus creating the conditions for a logarithmic integration over u. To explicitly see 
this, recall that we consider the weak-scattering regime, where the product of S'-matrices multiplying Af stl 
can be linearized. This allows us write SxuSuzSzy — SxuSuy — —Tuz — T^y + T^y — —2T(z). (Once again, 
the dominant contribution has been generated by the ‘real’ term.) We see that the net scattering amplitude 
in this approximation is independent of the intermediate dipole size u. Accordingly, when integrating over 
u, within the range limited by r and 2 , we find a STL, as anticipated. After also including the LO term and 
the NLO one enhanced by the DTL, one finds that the NLO BK equation in the strongly ordered region (H 
reduces to 


dr(r) 

dr 





( 5 ) 


It is now clear that, if the daughter dipoles are allowed to become sufficiently large, the NLO contributions 
enhanced by large transverse logarithms become comparable to, or larger than, the LO one. In that case, 
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the present perturbative expansion cannot be trusted anymore. To be more explicit, consider a single step 
AT in the evolution with the following, simple but physically meaningful, initial condition 


r(r) = 



for <C 1 

for r'^Ql » 1. 


The z-integration in Eq. ([^ becomes logarithmic and gives 
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(7) 


This shows that, for sufficiently small rQg, such that dg ln^(l/r^(5g) > 1, the NLO correction becomes larger 
than the LO term and the perturbation series is unreliable. In particular, the NLO correction is negative 
rendering the solution unstable, as indeed observed in numerical solutions [iSHSO] . 


3. Resumming the large collinear logarithms 


The large NLO corrections that we have singled out in the previous section are the lowest-order examples 
of collinearly-enhanced radiative corrections, which occur to all orders and spoil the convergence of the 
perturbation theory. When the separation between the transverse scales of the projectile and the target 
is large enough, as is actually the case in the DIS kinematics at HERA, the higher-order terms of this 
type become more important than the pure NLO terms (i.e. the contributions of 0{a1) which are not 
amplified by any transverse logarithms). From now on, we shall focus on this situation, discarding the pure 
dg NLO corrections, but focusing on the resummation of the large transverse logarithms to all orders. In 
this section, we consider both the single and double collinear logarithms, thus following and expanding our 
recent results in |50j . In the next section, we shall explain how the running coupling corrections can be 
included in this scheme. 

In Ref. [SD] we have devised a strategy for resumming double-logarithmic corrections to the BK equation 
to all orders. Our main observation was that these corrections are generated by the diagrams common 
to the BFKL and DGLAP evolutions — i.e. the Feynman graphs of light-cone perturbation theory in 
which the successive gluon emissions are strongly ordered in both longitudinal momenta and transverse 
momenta (or ‘dipole sizes’) — after enforcing the additional constraint that the emissions must also be 
ordered in lifetimes (or, equivalently, in light-cone energies |56)L Concerning the single collinear logarithms, 
it is intuitively clear that they must represent DGLAP-like corrections to BFKL, ‘small-a:’, emissions. For 
instance, the effect of order a^AYp^ visible in Eq. Q is the result of a sequence of two emissions: one 
small-a: emission (in the double logarithmic regime) yielding a contribution oc dgATp, and a DGLAP-like 
emission, characterized by strong ordering in dipole sizes (see eq. Eq. (H) and which gives an effect of 
order Here, p = is the large logarithm generated by the transverse phase-space. This 

scenario is corroborated by the following observation: the numerical coefficient Ai = 11/12 in front of the 
STL in Eq. ([^ can be recognized as the second-order term in the small lo expansion of the relevant linear 
combination of DGLAP anomalous dimensions: 


= f dz; 
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( 8 ) 


Recalling that one needs one factor of l/w in order to generate a small-x logarithm AT = ln(l/x), one sees 
that the NLO effect ~ d^AYp^ is indeed produced by combining the singular (l/w) piece of one emission 
with the first non-singular piece (Ai) of another one. This discussion also instructs us about the strategy 
to follow in order to resum such STLs to all orders: it suffices to include this piece Ai as an anomalous 
dimension, i.e. as an extra power-law suppression, in the evolution kernel previously obtained in Ref. [50]. 
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Figure 1: Running coupling for various schemes and configurations, (a) As a function of the daughter dipole 
size \x — zj, with </) = 0 the angle between the parent dipole x — y and the daughter one x — z. (b) The 
same with (j) = 7r/6. (c) As a function of the angle (p for fixed daughter dipole size |a; — 2 ;| = 1.5. Black 
(continuous): The minimal dipole scheme as defined in Eq. (11). Red (dashed): The “fac” scheme as given 
in Eq. (12). Blue (dotted): The Balitsky scheme [35]. In all cases the parent dipole size is \x — y\ = 1, the 
coupling is smoothly frozen at the value 0.7 and Aqcd = 0.2. 


We are thus led to the following, collinearly-improved, version of the BK equation, 
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where the overall kernel is written as a product of three factors: the familiar dipole kernel which appears 
already at leading order, the ‘DLA kernel’, resuming the double collinear logs to all orders |5D] 
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( 10 ) 


evaluated at p = \/lfczflfyzr, with = ln[(a; — 2 :)^/r^], and a new factor, which features the exponent 
idgAi (the positive sign in the exponent is taken when \x—y\ < min{|a: — 2 ;|, |y— 2 ;|} and the negative sign 
otherwise), which expresses the contribution of the single collinear logarithms. 

From the above discussion, is should also be clear that the present resummation of STLs is only partial: 
it refers to the particular class of such corrections which are generated by the first non-singular piece in the 
expansion in Eq. The higher terms in this w-expansion will produce single collinear logarithms too, but 
only starting at higher orders in perturbation theory (NNLO or higher). At the level of the BEKL equation, 
more complete resummations of the single logarithms have been devised in |52H54] , but so far it is not clear 
how to extend these resummation schemes to a non-linear evolution equation like BK. 

Returning to Eq. ([^ , the tilde symbol in f^y is intended to remind that this is truly a suitable analytic 
continuation of the dipole amplitude which coincides with the physical quantity T^y only for p < Y. Eor 
p > Y, the physical amplitude can be obtained by either solving an equation non-local in K, or by matching 
onto the solution to the DGLAP equation [SDj. However, explicit numerical studies at DLA level have shown 
that the solution f^y to Eq. © remains very close to the actual physical amplitude, including for p > Y. Eor 
this reason, we shall ignore this subtlety (and the related issue of the resummation in the initial condition) 
for the purpose of the fits to be constructed in the next section. We shall return to a more detailed study 
of these issues in a forthcoming publication [53] . 


4. Prescriptions for the running of the coupling 

The last source of potentially large NLO corrections to the BK equation are the running coupling 
corrections, i.e. the logarithmic terms proportional to b in the SI term in Eq. 0- Such terms can grow large 
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when the scales in their arguments are very disparate. More precisely, the first logarithm can be problematic 
when r is much smaller or much larger than l//i, while the second when the soft gluon at z is collinear to 
the quark or the antiquark composing the parent dipole. We need to choose /r in such a way to cancel these 
potentially large logarithms, which could otherwise spoil the convergence of the perturbative expansioiQ It 
is clear that there is not a unique choice, but in QCD one usually expects the hardest scale to determine 
the running of the coupling. Indeed, a quick inspection shows that the smallest dipole prescription 

dmin = ds(rmin) with Tmin = min{|a;-y|, |a;- 2 !|, |y- 2 !|} (11) 


cancels the large logarithms in all kinematic regions. 

Another possibility is to choose p so that all the terms with coefficient vanish. Given that in the 
current work we neglect all finite (i.e. not enhanced by a large logarithm) aj. terms, this looks like what is 
called the “fastest apparent convergence” (fac) scheme |60H62j . It is convenient in the sense that one is left 
with just the leading term in dg. We find that 


<^fac 


1 ^ [x-zf - {y-zf a,{\x-z\) - a,{\y-z\) 

as{\x-y\) {x-yY as{\x-z\)as{\y-z\) 


( 12 ) 


and it is an easy exercise to show that it reduces to the minimal dipole choice ds(r'niin) in all limits where 
one of the three dipoles is much smaller than the other two. 

In this work, we shall use both above schemes. Let us add that the most popular prescription, widely 
used so far in phenomenological applications, is the one due to Balitsky |32], and reads 


OBai = as{\x-y\) 


1 + 


as{\x-z\) - aY\y-z\) (y-s{\x-z\){y-zY - as{\y-z\){x -z) 
as{\x-z\)as{\y-z\) {x-yY 
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(13) 


but it will not be adopted here for a number of reasons. First, it is based on an extrapolation to all orders 
of a coordinate space kernel which includes the <5qb terms above as well as (YJY corrections. At this order, 
we would also expect corrections proportional to the two-loop beta function. Second, even though it also 
reduces to ds(7’min) in the extreme kinematical limits, it does that very slowly for large daughter dipoles (in 
certain configurations) and this leads to an unphysically small coupling in a large region of phase space, as 
can be seen in the respective plots in Fig. Finally, and perhaps as a result of the above drawbacks, when 
used in fitting the DIS data, it gives a much worse fit than Eqs. (11) and (12) and with fit parameters which 
take somewhat unnatural values. 


5. Fits to the HERA data 


We now turn to the description of the HERA reduced cross-section measurements using the resummed 
BK equation. To this aim several ingredients first have to be specified. 


Initial condition. We must fix the initial condition for the dipole amplitude at some Yoj which afterwards 
will be evolved towards higher rapidities using Eq. We consider two choices: the simple parametrisation 
of the Golec-Biernat and Wiisthoff (GBW) [5] type 


T{Yo,r) 


11 — exp 



(14) 


and the running-coupling version of the McLerran-Venugopalan (rcMV) model [58] 


T{Yo,r) 


11 — exp 



1 -1- In 



Pi 



(15) 


^It is rather important to point out here that p should cancel only these logarithms and not those discussed earlier which 
are of different physical origin. Of course one can proceed to such a choice and cancel all the NLO logarithms, but the result 
will be extremely unstable w.r.t. small variations of p. 






















init 

cdt. 

RG 

schm 

sing. 

logs 

per data point 

parameters 

^ red 

CC 

•^red 

Fl 

Rp[fm] 

Qo[GeV] 

a 

P 

GBW 

small 

yes 

1.135 

0.552 

0.596 

0.699 

0.428 

2.358 

2.802 

GBW 

fac 

yes 

1.262 

0.626 

0.602 

0.671 

0.460 

0.479 

1.148 

rcMV 

small 

yes 

1.126 

0.565 

0.592 

0.707 

0.633 

2.586 

0.807 

rcMV 

fac 

yes 

1.228 

0.647 

0.594 

0.677 

0.621 

0.504 

0.541 

GBW 

small 

no 

1.121 

0.597 

0.597 

0.716 

0.414 

6.428 

4.000 

GBW 

fac 

no 

1.164 

0.609 

0.594 

0.697 

0.429 

1.195 

4.000 

rcMV 

small 

no 

1.093 

0.539 

0.594 

0.718 

0.647 

7.012 

1.061 

rcMV 

fac 

no 

1.132 

0.550 

0.591 

0.699 

0.604 

1.295 

0.820 


Table 1: values of the fitted parameters entering the description of the HERA data. The fit includes 

the 252 (Tred data points. The quoted for and Fl are obtained a posteriori. 


It is worth noticing that, as dictated by collinear physics, there is no anomalous dimension in the above 
initial conditions. The extra parameter p determines the shape of the amplitude close to saturation and its 
approach towards unitarity. 


Running coupling. We consider the two prescriptions given by Eqs. (11) and (12). Eor the explicit expression 
of the strong coupling in coordinate space in terms of r we introduce a fudge factor as in [22) . namely 


o-sir) = 


1 


bNt In [4C2/(r2A^^)] ’ 


(16) 


with h]si^ = (llA^c — 2A^f)/127r. This fudge factor is also included in the rcMV type initial condition in (15). 
The A^f-dependent Landau pole is obtained by imposing as(Mj) = 0.1185 at the scale of the Z mass m 
and continuity of ttg at the flavour thresholds, using rUc = 1.3 GeV and = 4.5 GeV. To regularise the 
infrared behaviour, we have decided to freeze at a value agat = 1 and we have checked explicitly that 
reducing this down to, for example, 0.7 does not affect the fit in any significant manner. 

Note that we do not include any form of resummation or matching for Inl/r^ > E, as introduced in 
m, in these initial conditions. One of the reasons for not doing so is that the extra factor in the initial 
condition can always be reabsorbed in a re-parametrisation. Eurthermore, a proper matching at small dipole 
sizes, suited for phenomenological studies, would require a careful treatment of the small-dipole region. In 
that respect, the resummed BK evolution is expected to perform a better job than a fixed matching with 
a fixed asymptotic behaviour. We leave a better treatment, e.g. a genuine matching to DGLAP evolution, 
for future work. 

Rapidity evolution. Of course this is determined by the resummed BK equation given in ([^ . Here, we again 
consider two separate cases, one in which the evolution resums only the leading double logarithms and one 
in which it also includes the single ones. 

From the dipole amplitude to observables. Once we have the dipole amplitude for all rapidities and dipole 
sizes, we use the standard dipole formalism to obtain the physical observables: 


= 2TrRlY^ J d^r j dz |4'[f^(r, z; ^(In 1/%, r). 


where the transverse and longitudinal virtual photon wavefunctions read 
,(/)/'», r12^|2 _ „2'^emA^c ! 


\d>^^\r,z-QX = <^^dQ^z\l-zYKt,{rQf) 


2tF- 


{ V' + (1 - Q}Kl{rQf) + m}Kl{rQf)] . 
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(17) 

(18) 
(19) 




















init 

cdt. 

RG 

schm 

sing. 

logs 

xVnpts for 

50 

100 

200 

400 

GBW 

small 

yes 

1.135 

1.172 

1.355 

1.537 

GBW 

fac 

yes 

1.262 

1.360 

1.654 

1.899 

rcMV 

small 

yes 

1.126 

1.170 

1.182 

1.197 

rcMV 

fac 

yes 

1.228 

1.304 

1.377 

1.421 

GBW 

small 

no 

1.121 

1.131 

1.317 

1.487 

GBW 

fac 

no 

1.164 

1.203 

1.421 

1.622 

rcMV 

small 

no 

1.093 

1.116 

1.106 

1.109 

rcMV 

fac 

no 

1.131 

1.181 

1.171 

1.171 


Table 2: Evolution of the fit quality when including data at larger (in GeV^). 


In the above we have introduced the customary notation = z(l — z)Q’^ + m^, i/ = x{l + and 

we have assumed a uniform distribution over a disk of radius in impact parameter space. The sum in ( |17[ ) 
runs over all quark flavours and we will include the contributions from light quarks with mu^d,s = 100 MeV 
as well as from the charm quark with ijif. = 1.3 GeV. From the longitudinal and transverse cross-sections, 
we can deduce the reduced cross-section and the longitudinal structure function as 


<^red 




dTr^Oem 

dTr^Oem 


2(1-1/) 7*P 

l + (l-y)2'"L 


( 20 ) 

( 21 ) 


When the quark masses, the value of the strong coupling at the Z mass and its frozen value in the 
infrared have been fixed, we are left with 4 free parameters according to our choice of initial condition: Rp 
the “proton radius”, Qq the scale separating the dilute and dense regimes, Ca the fudge factor in the running 
coupling in coordinate space, and p which controls the approach to saturation in the initial condition. 

We have htted these parameters to the combined HERA measurements of the reduced photon-proton 
cross-section [38) . Since the BK equation is applicable only at small-cc, we have limited ourselves to the 
region x < 0.01. We note that since Eq. 0 probes dipoles at the rapidity In l/i/, the exact cut we impose 
is Xc < 0.01 since the most constraining cut comes from the charm, the most massive quark we include in 
our model. Accordingly, our initial condition for the BK evolution corresponds to i = 0.01. Furthermore, 
since we do not expect the BK equation to capture the full collinear physics, we impose the upper bound 

< Qmax- Hy default we will use Q^ax = 50 GeV^ but we will also give results for extensions to larger 

. In the default case we have a total of 252 points included in the fit. We have added the statistical and 
systematic uncertainties in quadrature]^ 

The results of our fits for the 2^ = 8 cases, depending on the initial condition, the running coupling 
prescription and the inclusion or not of single logarithms in the kernel, are presented in Table. The table 
includes the parameter values obtained from fitting the (Jred data and, besides the fit it also indicates 
the obtained a posteriori for the latest and Fl [SS] measurements. These results deserve a few 

important comments. 

(i) In general, the overall quality of the fit is very good, reaching per point around 1.1-1.2. 

(ii) Apart from a few small exceptions (see below), all the parameters take acceptable values of order one. 
Note that we have manually bounded p between 0.25 and 4. Whenever it reached the upper limit, 
larger values only led to minor improvements in the quality of the fit. 


3A more involved treatment of the correlated systematic uncertainties leads to similar results with slightly worse P®*' 
points (about 0.04). 
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Figure 2: Description of the HERA data obtained by the fits using the rcMV initial condition. Each box 
corresponds to a given value of as indicated (in GeV^) in the top-right corner. For each fit we plot 
the ratio of the prediction to the central experimental value. The (green) band represents the experimental 
uncertainty. 
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(iii) 


(iv) 


(v) 


(vi) 


The two initial conditions give similar results, with a slight advantage for the rcMV option, which is 
likely due to the extra parameter. Note that for a standard MV-type of initial condition T{YQ,r) = 
{1 - exp[-(r^(5Q/4 [c + ln(l + l/rA)])^]}^/^, we have not been able to obtain a psr point below 1.3 
and the parameters, typically c or p, tend to take unnatural values. 

As far as the running-coupling prescription is concerned, the smallest dipole prescription given in 
Eq. @ tends to give somewhat better fits than the “fac” prescription given in Eq. ( [T^ . This can be 
seen as an estimate of subleading corrections (including the pure a'^ NLO terms) that we neglect in 
the present fit. Note also that we have not been able to reach a fit of equivalent quality and robustness 
with the Balitsky prescription. 

The resummation of the single logarithms tends to yield slightly larger values for , but the difference 
is too small to be significative (at least, without performing a full NLO analysis). Perhaps more 
significantly, this resummation leads to more physical values for some of the parameters, especially Cq 
for the smallest dipole prescription and p for the GBW initial condition. These findings are consistent 
with the expectation that, once properly resummed, single logarithms should have only a modest 
impact. Recall however that their resummation is a crucial step towards a full NLO fit — failing to 
do so could lead to instabilities similar to those observed when double logarithms are not resummed. 
The fit remains stable when varying the parameters we have imposed by hand. For example, using 
Osat = 0.7 instead of 1 has no significant effect on the fit. Varying the light quark masses within the 
rather wide range 0 < mu,d,s < 140 MeV only slightly changes the quality of the fit. For instance, 
taking one of our best fits (rcMV initial condition, the smallest dipole prescription for the running of the 
coupling, and resummation of the single logarithms), we have found x^ = {1-180, 1.153, 1.126, 1.159} 
when choosing m„ ^ = {0, 50, 100, 140} MeV, respectively. This lack of sensitivity to the light quark 
masses is likely a consequence of saturation, which reduces the dependence of the DIS cross-section to 
very large dipole fluctuations. (The corresponding amplitudes reach the unitarity, or ‘black disk’, limit 
r = 1, so they are independent of the size r of the dipoles fluctuations, as regulated at low Q'^ by the 
quark masses; see also the discussion of Fig. below.) Also, we have obtained an equally good fit with 
the slightly larger value m-c = 1.4 GeV for the mass of the charm quark, although the quality started 
deteriorating for significantly larger values rric >1.6 GeV. Very similar findings have been reported 
for the saturation fits in [50]. 

Trying to extend the fit to larger shows an interesting behaviour as seen from Table While the 
obtained using the GBW initial condition increase when including higher-Q^ data, the fits using 
the rcMV initial condition remain stable. We suspect that this is due to the fact that this choice of 
initial condition stays closer to the expected physics at high 

In Fig. 1^ one can see the quality of our fit and the extracted values of the evolution parameter As = 
dlnQg/dV. In Fig. we show the value of the saturation momentum in the (a;, (5^)-plane on top of the 
data points as well as a few selected initial conditions for the fit. Note that amplitudes which a priori 


(vii) 


have different functional forms, cf. Eqs. (14) and (15), look nevertheless quite similar in shape (at least in 
double-logarithmic scale) when plotted for the particular values of the parameters that are selected by the 
fits. 

To conclude, this work can be seen as the first description of small-a; DIS data which includes mandatory 
perturbative QGD ingredients in that region: leading-order small-a; evolution, the resummation of large 
transverse logarithms, and saturation correction^ The dipole amplitude obtained from our fits to inclusive 
DIS can in principle be used to compute several other observables, like particle multiplicity in hadronic 
collisions, the diffractive structure functions, the elastic production of vector mesons, or the forward particle 
production in heavy-ion collisions. This is certainly not the end of the story: beyond what we have included 
here, it would be interesting to add the pure aj. NLO corrections to the BK evolution kernel, thus obtaining 
a genuine resummed NLO-BK fit, and to perform a proper matching between this small-cc evolution and a 
DGLAP-like evolution at large and large x. These steps go beyond the scope of the present paper and 
are left for future studies. 


^It would be an interesting exercise to see what happens if one switches off the non-linear corrections in Eq. |^. Given our 
asymmetric choice of frame, justified by saturation physics, this may however require extra work. See also ESj. 
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Figure 3: Left: value of the saturation momentum, defined for each rapidity as 2/rs{Y) with T(rs(F), y) = 
1/2. For comparison, we have overlaid the experimental data points from the HERA dataset. Right: plot 
of the corresponding initial conditions for the rapidity evolution. 
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